Interference between competing pathways in atomic harmonic generation 
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We investigate the influence of the autoionizing 3s3p 6 n£ resonances on the fifth harmonic gener- 
ated by 200 - 240 nm laser fields interacting with Ar. To determine the influence of a multi-electron 
response we have developed capability within time dependent 7?-matrix theory to determine the 
harmonic spectra generated. The fifth harmonic is affected by interference between the response of 
a 3s electron and the response of a 3p electron, as demonstrated by the asymmetric profiles in the 
harmonic yields as functions of wavelength. 

PACS numbers: 32.80.Rm, 31.15.A-, 42.65.Ky 



Harmonic generation (HG) is one of the fundamen- 
tal processes that can occur in light-matter interac- 
tion. The process is of strong technological interest, 
since it lies at the heart of ultra-short light pulse gen- 
eration [l|. However, HG can also be used as a mea- 
surement tool, for example, to study molecular changes 
on the sub- femtosecond timescale 0]. Recent articles 
have postulated various other applications of high har- 
monics, including electron interferometry and the tomo- 
graphic imaging of electron orbitals [3] . HG has also been 
used experimentally to demonstrate the importance of 
electron-correlation in ultra-fast molecular dynamics Q 
and multi-electron dynamics in atomic systems 

In order to probe the influence of multi-electron dy- 
namics on HG theoretically, we find it more expedient to 
consider atomic systems. Although detailed calculations 
for molecules are possible using sophisticated codes, it 
is much more difficult to identify the multi-electron dy- 
namics within a molecular system than within an atomic 
system. Even so, many methods developed for the deter- 
mination of harmonic spectra in atoms rely on the single 
active electron approximation [(| 0|- Several methods 
have been developed for the description of two-electron 
systems , but these methods are not easy to extend 
to general atoms. New methods are thus needed to gain 
insight into the influence of multi-electron dynamics on 
harmonic spectra in general multi-electron atoms. 

At Queen's University Belfast we have developed time- 
dependent iZ-matrix (TDRM) theory over the last five 
years expressly to study the behavior of multi-electron 
atoms in short, intense light fields [10|. The advantage 
of this method in studying multi-electron dynamics was 
shown in a pump-probe study of C + ll]: by varying the 
delay between the pump and probe pulses, oscillations in 
the total ionization yield were observed, which could be 
linked to the correlated motion of two electrons. Because 
of the interest in using harmonic spectra to observe mult i- 
electron dynamics, we have developed capability within 
the time-dependent iZ-matrix approach to determine har- 
monic spectra for a single atom. 

As a first demonstration of this capability, we report 
on the influence of the 3s3p 6 n£ Rydberg series on single 
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FIG. 1: (Color online) The harmonic yield in Argon is de- 
termined by two competing pathways: absorption of at least 
three photons can excite a 3p electron (red solid line) into the 
continuum, while absorption of at least five photons can ex- 
cite a 3s (blue dashed line) electron to the np Rydberg series 
converging onto the 3s3p 6 threshold. 



atom harmonic spectra in Ar. In particular, we focus on 
the wavelength range for which the fifth harmonic coin- 
cides with these resonances: 200 - 240 nm. Recent stud- 
ies have focused on the parallels between photoionization 
and HG [IH, [l3|, where, in particular, the Cooper min- 
imum in Ar was investigated. However, minima in the 
photoionization cross section can also occur near window 
resonances, where resonances actually lead to a decrease 
in the photoionization cross section. The 3s3p 6 n£ res- 
onances are examples of window resonances [14| . Their 
decay towards 3s 2 3p 5 e£' continua requires a good descrip- 
tion of dielectronic repulsion. The harmonic yield will be 
determined by two interfering pathways, (Fig. [1}: exci- 
tation of a 3s electron to np, and excitation of a 3p elec- 
tron into the continuum. A multi-electron calculation is 
essential to describe the harmonic yield accurately. This 
study thus provides a good test case. 

To describe the extension to compute a system's har- 
monic response, we first give a brief overview of TDRM 
theory. A more thorough discussion can be found in 
The time-dependent Schrodinger equation for an atom 
containing (N + 1) electrons, and nuclear charge Z, is 



d 

i— ¥(X w+1 ,t) 



H{t)*(X N+1 ,t). 



(1) 



2 



The Hamiltonian, H , contains both the non-relativistic 
Hamiltonian of the N + 1-electron atom or ion in the ab- 
sence of the laser field and the laser interaction term. The 
laser field is described using the dipole approximation in 
the length form, and is assumed to be linearly polarized 
and spatially homogeneous. This form provides the most 
reliable ionization yields when only a limited amount of 
atomic structure is included fl5| . 

To propagate the wave function in time, time is dis- 
cretized into steps of At such that 

t q = qAt, g = 0,l,2,.... (2) 

Then, using the unitary Cayley form of the time evolution 
operator, exp [— itH (t)], the wavefunction ^ in Eq. (TTJ) 
at tq+i can be written in terms of the wavefunction at t q 

[H (t q+1/2 ) -E] ¥(X w+ i,i,+i) =®(X N+1 ,t q ), (3) 

where 

&(X N+1 ,t q ) = - [H(t q+1/2 )+E]*(X N+1 ,t q ), (4) 

and H (^+1/2) is the time-dependent Hamiltonian at the 
midpoint. E is given by 2i/At. In equations ([3]) and (|4j) 
Xn+i = ricri,^^, ••• , rjv+iCTjv+i where r i a l are the 
space and spin coordinates of the i th electron. If the 
wavefunction at time t = t q is known, then Eq. (|3j allows 
us to propagate the solution forward in time. 

In i?-matrix theory, configuration space is partitioned 
into two regions: an inner and an outer region. In the 
inner region, all electrons are within a distance a m of 
the nucleus, and all interactions between all electrons are 
taken into account. In the outer region, one electron has 
moved beyond a distance a- m , and exchange interactions 
between this electron and the electrons remaining close 
to the nucleus can be neglected. The outer electron thus 
experiences only the long-range multipole potential of the 
residual ./V-electron core and the laser field. 

In the outer region, a ln < r^+i < a ou t, we have 

*(Xjv+i,f g+ i) = 

(Xat; iN+KJN+^r^Fp (rjv+i) , (5) 

P =i 

where <E>^ are channel functions which couple possible 
final states of the Ar + ion with spin and angular mo- 
mentum of the ejected electron. The F p are the radial 
wavefunctions which describe the motion of the ejected 
electron in the p th channel, where p is a channel index. 

We can now follow i?-matrix procedures [l6[ to evalu- 
ate Eq. ([3]) at the boundary a- ln as a matrix equation 

F (a in ) = R(a in )F (a in ) + T (a in ) , (6) 

in which the wavefunction at the boundary is described in 
terms of its derivative, F, plus an inhomogeneous vector, 
T, arising from the O term in Eq. ([3]). Equation © 
connects the inner and outer region wavefunction at the 
boundary a; n . 



Given an initial inner region wavefunction, R and T 
are evaluated at the boundary. Subsequently, they are 
propagated outwards in space up to a boundary, a ou t 
where it can be assumed that the wavefunction F has 
vanished. The wavefunction vector F is set to zero and 
then propagated inwards to the inner region boundary. 
Once the F- vector has been determined at each boundary 
point, the full wavefunction can be extracted from the 
i?-matrix equations. Equation ([3]) then allows us to step 
this solution forward in time. 

HG arises from the oscillating dipole moment induced 
by the laser field. This is given by the expectation value 
of the electric dipole operator, 

d (t) = <*(*) I -eZ|*(t)>, (7) 

where e is the electronic charge and Z the total position 
operator along the laser polarization axis. The power 
spectrum of the emitted radiation is then given by the 
Fourier transform of either the dipole moment, its veloc- 
ity or acceleration [17] . Many techniques involving high 
order harmonics use the acceleration representation as 
it gives better resolution for the high-order peaks fl8j . 
However, this has predominantly been verified for single- 
active-electron systems, whereas in the multi-electron 
TDRM method the length form of the dipole operator 
has proven to be far more reliable (lfjj . 

To describe Ar, we use the i?-matrix basis developed 
for single photon-ionization [13 ] with an inner region 
radius of 20 a.u. The set of continuum orbitals con- 
tains 60 _B-splines for each angular momentum, I, of the 
continuum electron. We use both the 3s 2 3p 5 2 P° and 
3s3p 6 1 P° states of Ar+ as target states. The descrip- 
tion of Ar includes all 3s 2 3p 5 e£ and 3s3p 6 e£ channels up 
to L max = 15. In the outer region we propagate R and T 
out to a radial distance of 1600 a.u. to prevent any reflec- 
tions of the wavefunction from the outer region boundary. 
The outer region is divided into sectors of width 2 a.u. 
and each sector contains 35 S-splines of order k = 9 per 
channel. We use similar pulse shapes to those used in 
[icj |: a 3 cycle sin 2 turn on/off for the electric field. Peak 
intensity is set at 2 x 10 13 Wcm~ 2 with varying numbers 
of cycles at peak intensity; the time step is 0.1 a.u. 

Figure [2] shows a typical power spectrum obtained for 
an incident pulse with a photon energy of 5.4 eV and 25 
cycles at peak intensity. As well as the three harmonic 
peaks, several other features are visible. The sideband 
structure of the harmonic peaks is due to the turn on/off 
of the laser pulse. In the photon energy region between 
12.1-15.2 eV, several peaks are observed arising from the 
3s 2 3p 5 n£ Rydberg series: 4s at 12.1 eV, the 5s/3d man- 
ifold at 14.5 eV and the 4d/6s at 15.1 eV. Higher man- 
ifolds are not separated at this pulse length. The third 
harmonic peak is clearly visible; it lies just above the 
3s 2 3p 5 threshold of Ar + . Its proximity to the threshold 
may mean that it still overlaps the Rydberg series lead- 
ing up to this threshold. The peak at 18.8 eV cannot 
be associated directly with any atomic structure but this 
energy corresponds to the energy of 3s 2 3p 5 4p plus one 
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FIG. 2: (Color online) The power spectrum produced by a 
25 cycle, 2 x 10 13 Wcm -2 laser pulse of photon energy 5.4eV. 
The features labelled are (a) First harmonic peak, (b) Effects 
due to 3s 2 3p 5 n£ Rydberg series, (c) Third harmonic peak, 
(d) Peak due to 3s 2 3p 5 4p states plus a photon, and (e) Fifth 
harmonic peak. Inset: The fifth harmonic region of the power 
spectra for two 25 cycle pulses of photon energy 5.36 eV (blue, 
dashed line), and 5.42 eV (red, solid line). The change in the 
profile of the peaks is typical of the presence of a resonance 
in between these two frequencies. 



5.4 eV photon. Finally, the fifth harmonic is visible at 
an energy of 26.9 eV. 

The power spectrum was compared with calculations 
carried out using the J?-matrix Floquet (RMF) method 
20]. Good agreement was found at 248nm for the 
first harmonic rate at intensities between 6.5 and 8.5 x 
10 12 Wcm -2 . The RMF calculations also obtained third 
and fifth harmonic rates within the 3d/5s manifold. In 
RMF theory these states are separate, but this requires 
unfeasibly long pulses in TDRM theory. Thus no com- 
parison is possible beyond the first harmonic. 

The main interest in the present manuscript is the fifth 
harmonic, which is strongly influenced by the 3s3p e n£ 
Rydberg series. Photoionization experiments and calcu- 
lations have demonstrated that the Rydberg series ap- 
pears as window resonances 00]. Figure [2] also shows 
part of the power spectra at two distinct photo ener- 
gies, 5.36 eV and 5.42 eV. The resonance profile in the 
harmonic spectrum changes significantly. This is typi- 
cal when a resonance lies in between these two photon 
energies. In order to assess these changes, the fifth har- 
monic rate is determined over a range of photon energies. 
For each photon energy, we integrate the harmonic rate 
across the bandwidth of the harmonic peak. This ap- 
proach yields the harmonic rates shown in Fig. [3] We 
can see resonances occurring at energies corresponding to 
the 3s3p 6 4p and 3s3p 6 5p states. Additional resonances 
may be present at higher photon energies, but the energy 
uncertainty makes their assessment difficult. 

In order to obtain resonance parameters, calculations 
are also performed with a reduced basis: a 3s 2 3p 5 single- 
target-state description of Ar + , and an Ar basis which 



FIG. 3: (Color online) The fifth harmonic rate for Ar in a 25 
cycle pulse of varying frequency, for a target with the full con- 
figuration (blue, solid line) and one without any contribution 
from the 3s3p l 'n/e£ series (red, dashed line). The resonances 
in the fifth harmonic rate correspond to the 3s3p 6 4p (26.9eV) 
and 3s3p 6 5p (28.3 eV) resonances. 



contains no 3s3p 6 n/s£ functions. This basis produces a 
smooth spectrum with a rate which can act as a back- 
ground. Subtraction of this background rate then re- 
veals a resonance profile due to the 3s3p 6 n£ Rydberg 
series, from which resonance parameters can be estab- 
lished. The resonance profile in Fig. @]was obtained for 
25 cycle pulses at an intensity of 2 x 10 13 Wcm -2 , giv- 
ing, within an accuracy of 10%, an asymmetry parameter 
q = -1.84 and line width, T = 0.27 eV. 

The resonance parameters differ significantly from 
those obtained in previous photoionization calculations 
[l9l |. The present resonance width is about a factor 4 
larger than obtained previously, 0.068 eV. This is not 
too surprising, since the present calculations are carried 
out only over a finite laser pulse, which gives rise to a 
broadening of the resonance. The effects of this broad- 
ening can be estimated by carrying out calculations for 
different pulse lengths. Table U indeed shows that the res- 
onance width decreases with increasing laser pulse length. 
In addition, the resonance may also broaden due to the 
intensity of the laser pulse. At higher intensities states 
can photoionize as well as autoionize. This leads to a 
reduction in the lifetime of the state, and hence a broad- 
ening. This effect is shown in Table [Til where resonance 
parameters are compared as a function of peak intensity. 

Figure [4] shows that at the fifth-harmonic stage, the 
resonance has a strong effect on the harmonic yields. 
The strong interference between the continuum and the 
resonance leads to an asymmetric profile with harmonic 
yields increased below and decreased above the resonance 
energy. Both the continuum and the resonant pathway 
to HG must therefore be accounted for. As these path- 
ways involve the response of different electrons, it shows 
that a multi-electron treatment of HG is necessary for Ar 
in the 200-240nm wavelength range. 

Table [J shows that the asymmetry parameter is not af- 
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FIG. 4: (Color online) The fifth harmonic resonance profile for 
Ar in a 25 cycle laser pulse of intensity 2 x 10 13 Wcm -2 . The 
data (blue points) is fitted by a Fano profile (red line) with 
asymmetry parameter, q = —1.84 and line width V = 0.27 
eV, and a linear background term. 

fected significantly by the pulse duration. On the other 
hand, Table HI1 shows that this parameter changes notice- 
ably with the intensity of the laser pulse. This knowl- 
edge is important for modelling the harmonic response of 
a macroscopic system and the appearance of resonances 
over a particular energy distribution. However, it is dif- 
ficult to analyze the parameter changes in detail. At an 
intensity of 5 x 10 13 Wcm -2 the ponderomotive shift of 
the ground state means that at the fifth-harmonic reso- 
nance, the third harmonic lies very close to the 3s 2 3p 5 
threshold, and this may affect the profile. 

In conclusion we have applied TDRM theory to study 
the appearance of the 3s3p e n£ autoionizing resonances in 
the fifth-harmonic for Ar. These window resonances show 
significant asymmetry within the fifth-harmonic spec- 
trum, indicating interference between pathways involving 
the response of different electrons. This study provides 
a first demonstration of the new capability within time- 
dependent R-matrix theory to obtain harmonic spectra 
in multi-electron atoms with electron-electron correlation 
fully included. This extension may enable further evalu- 
ation of whether and how evidence of multi-electron dy- 
namics can be found within harmonic spectra of atoms. 



In the future, we hope to extend the approach to longer 
wavelengths and higher intensities. This will, however, 
require an assessment of the different gauges that can 
be used to obtain the harmonic spectrum, and extensive 
code development. The TDRM approach has proven its 
merit for strong-field ioniztion at wavelengths of 390nm 
at 4 x 10 14 Wcm~ 2 [l6[. It will be of interest to export the 
determination of harmonic yields to the newly developed 
i?-matrix with time-dependence codes (2lj j . which may 
be better suited for longer wavelengths. 
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TABLE I: The asymmetry parameters, line widths and posi- 
tions for the resonance profiles due to 2 x 10 13 Wcm -2 laser 
pulses of various lengths accurate to within ~10% 



Peak cycles q parameter Line width (eV) Position (eV) 



10 


-1.83 


0.50 


26.99 


15 


-1.82 


0.40 


27.01 


20 


-1.83 


0.30 


27.00 


25 


-1.84 


0.27 


27.01 



TABLE II: The asymmetry parameters, line widths and po- 
sitions for the resonance profiles due to 25 cycle laser pulses 
of various intensities, accurate to within ~10% 



Intensity q parameter Line width Position 



Wcm" 2 ) 




(eV) 


(eV) 


2.0 


-1.84 


0.27 


27.01 


3.0 


-1.79 


0.30 


27.08 


5.0 


-1.31 


0.32 


27.27 
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